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We show that the dipole-dipole interatomic forces induced 
by a single off-resonant running laser beam can lead to a self- 
bound pencil-shaped Bose condensate, even if the laser beam 
is a plane- wave. For an appropriate laser intensity the ground 
state has a quasi-one dimensional density modulation — a 
Bose "supersolid". 
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Electrostriction is the tendency of matter to become 
compressed in the presence of an electric field [Q. In 
an optically trapped atomic Bose-Einstcin condensate 
(BEG) it is provided by the gradient of the incident elec- 
tric field (one-body dipole forces). However, it can occur 
even if the external fields are homogeneous: the elec- 
trostriction effects of dipole-dipole interatomic forces in- 
duced by certain configurations of far-off-resonant laser 
beams are capable, for sufficiently high laser intensity, of 
self-binding a BEG within a wavelength of the laser (in 
the near-zone) |^]. Thereby, a physical situation can be 
realized which has analogies to self-gravity . 

Here we show that even in the simplest case of a single 
homogeneous, circularly-polarized, far-off-resonant laser 
beam a BEG can be self-bound by electrostriction, due 
to the retarded part of the dipole-dipole interaction. 
The self-bound condensate will be pencil shaped, aligned 
along the direction of propagation of light, in principle 
even for arbitrarily small intensity. As the laser intensity 
increases, the electrostriction can give rise to a remark- 
able density modulation of the quasi-one dimensional 
condensate ground state. Such a condensate bears simi- 
larity to a "supersolid", i.e., a long-range crystalline- like 
density modulation imposed upon a superffuid by inter- 
particle forces The formation of such structures is 
associated with a strong enhancement of the elastically 
scattered field (and, consequently, a large radiation pres- 
sure) and with suppression of the heating due to spon- 
taneous Rayleigh scattering. The relation of our effects 
to (stimulated) Brillouin scattering, collective atomic re- 
coil and superradiant Rayleigh scattering is briefly 
discussed. 

We consider the case of a far-off-resonant circularly 
polarized laser beam propagating along the positive z 
direction with wavevector q (Fig. 1 — inset). The dipole- 
dipole induced interatomic potential |q] then becomes 
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Here x, y, z are the components of the interatomic sep- 



aration r and hTray 



^q^I/GiTceQ is the single-atom 



(spontaneous) Rayleigh scattering rate, which is propor- 
tional to the laser intensity / and to a^, i.e. to the satu- 
ration parameter that scales as the inverse square of the 
detuning from the nearest atomic resonance |3j. We as- 
sume in the following that the condensate contains many 
atoms per cubic wavelength to ensure the validity of the 
mean-field approach. The condensate is taken to be at 
zero temperature and the rate of heating due to incoher- 
ent Rayleigh scattering is shown to be negligible. 

The Gross-Pitaevskii equation [|l0| for the order pa- 
rameter 5'(r,f) can be obtained starting from 
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Here the total mean-field energy functional i?tot = 
i?kin + H^o + Hdd + Hscat consists of: (a) the ki- 
netic energy, Hmu = J dr(S^/2TO)| V*p; (b) the 
harmonic-trap energy, that can be due to the focus 
size of the beam, iJho ~ J drVhon{r), where n ~ 
|\E'(r,t)|^ is the atomic density and Vho — + 
muj'^z'^/2 is the harmonic cylindrically symmetric po- 
tential, where = + y'^; (c) the mean- field energy 
due to the short-range (r~^) van der Waals interaction, 
which will be treated within the delta-function pseudo- 
potential approximation, iJscat = {2'Kafi'^ /m) J(irn(r)^, 
where a is the s-wave scattering length, and (d) the 
electromagnetically-induced mean-field energy, i/^d = 
(1/2) J drdr'Vddi'T — r')n(r)n(r') which corresponds to 
the electrostriction energy of the medium Q to lowest 
order in the an/eg. The mean radii of the condensate 
can be estimated using the Gaussian ansatz \E'(p, z) — 

exp (-pV2t«r - ^V2w2) /n^/^wV^Wr. The varia- 
tional parameters Wz and Wr, which are related to the 
mean quadratic radii through (Ax)^ = (x^) = and 
(Az)^ = (z^) — wl/2, are obtained by minimizing TJtot- 
A self-bound condensate is then found to exist in the 
Thomas- Fermi (TF) limit (negligible kinetic energy) . Ex- 
ternal confinement is not necessary and the values of Wr 
and Wz are finite due only to the presence of a plane-wave 
laser. The laser intensities are below the threshold of the 
instability caused by the 1 /r^ part of the dipole-dipole 
potential (|l|) (static field limit) / < 12TrTi'^ ce^a / ma'^ , 
and can be, in principle, arbitrarily small. Here and 
henceforth we normalize the intensity to Snh^ ce^a / ma'^ 
and the spatial scales are normalized to the laser wave- 
length Xl — 271 /q. 
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Figure 1 shows the variational values of the mean con- 
densate radii in the TF-limit for a plane wave laser and 
in the absence of external confinement. For / ^ 1 the 
radii of the condensate are given by Aa; — 0.1125/^^/^ 
and Az — 0.7847 /~^. We note that for laser intensity 
/ > 0.1 the condensate is strongly confined in the radial 
direction with Ax ~ 0.2, and less confined in the longi- 
tudinal direction, with typical size larger or of the order 
of the wavelength, Az > 0.6. 
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FIG. 1. Inset — the laser beam and condensate geome- 
try. Mean radii (normalized to Xl) of the self-bound con- 
densate versus the intensity / (normalized to SnfV^ceQa/ma^) 
in the TF-limit. Long-dashed (dotted) lines represent the 
expectation value of Az (Ax) obtained from the Gaussian 
ansatz. Circles (diamonds) represent the expectation value of 
Az (Ax) calculated from the quasi- ID ansatz (3). 

The large extension of the condensate along the lon- 
gitudinal (i) axis means that the oscillatory long-range 
behavior of the potential (|^) can become manifest in this 
direction. We study the effect of these oscillations using 
the ansatz 



*(p, z) = exp (-pV2w^) /vr^/^w^ 
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in order to reduce the problem to a quasi-one dimen- 
sional (ID) one. We then obtain ip{z), as well as Wr, 
by minimizing numerically the mean- field energy TJtot- 
The numerical solutions for the ground state, based on 
(^, yield radial and longitudinal sizes of the self-confined 
condensate (i?ho = 0) that are in good agreement with 
those based on the Gaussian ansatz, as displayed in Fig 
1. 

The electromagnetically-induced mean-field energy, 
expressed through the linear density niD = takes 
the ID-reduced form H^^ = (1/2) / dzdz'Vj^^{z - z') 
niT>{z)niY){z'), where the induced dipole-dipole poten- 
tial, reduced to ID, is given by 
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The ID-reduced interaction contains an attractive singu- 
lar part Kid°(^)L,„g. = -[KCray/q^'wl) 5{z). For / < 1.5 
the singular part of the reduced interaction is balanced by 
a similar repulsive term V^^^^{z) = {2ah^ /'m'w'f)5{z) pro- 
portional to the (positive) s-wave scattering length. Fig- 
ure 2 shows the non singular part of V^^^P for three differ- 
ent values of the variational parameter Wr ■ For large val- 
ues of z the ID-reduced potential oscillates as V^]^{z) k, 
— {3hTray/2q) cos(gz)^/z, giving rise to characteristic at- 
tractive logarithmic singularities in the Fourier transform 
at fc^ = and = ±2q (Fig. 2b). 
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FIG. 2. (a) the induced ID dipole-dipole interaction in 
units of hTray, as a function of z (normalized to Xl), (b) 
its Fourier transform. The dashed, solid and dotted lines 
correspond to increasing values of the variational parameter 
Wr — 0.1, 0.2 and 0.3, respectively. 

In the TF-limit, the Bogoliubov dispersion relation for 
the condensate |l^ is proportional to the square root of 
the Fourier transform of the interparticle potential. In 
a sense then, the = ±2q singularities give rise to a 
roton-like minimum (cf. He-II) in the dispersion relation. 
When the roton minimum touches the zero frequency axis 
there is no energy cost to the system adopting a density 
modulation in its ground state, and a 'supersolid' is cre- 
ated § (see below). 

A new feature appearing in the TF-limit is the for- 
mation of single condensate droplets of size less than 
~ .25Al located at an approximately regular distance 
of Al/2. The variational values for ground state energies 
based on (^ are substantially below the ones obtained by 
the Gaussian ansatz, as displayed in Fig. 3a, and indicate 
that such density modulation is quite stable as its energy 
reduction is of order of the Gaussian variational energy. 
The effect of density modulation arises as the field back- 
scattered by the condensate interferes with the incident 
field, thus modulating the total laser intensity and conse- 
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quently creating a series of attractive dipole traps along 
the longitudinal axis that are separated by half a wave- 
length (see thick curve in Fig. 4). In this limit, the single 
condensate droplets become more and more confined and 
their number is reduced, when approaching the threshold 
of the l/r^ instability (/ = 1.5). 
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FIG. 3. (a) TF mean-field energies per particle 

(normalized to tiEr) versus the intensity / (normalized 
to cel^a / mo?) as obtained from the Gaussian ansatz 
(solid-line) and from the ansatz (3) (diamonds), (b) Kinetic 
energy effects: longitudinal equilibrium densities for different 
values of r) and for 7 = 1. A radial external confinement is 
used to keep the radial size of the condensate constant and 
equal to that obtained for -q very large without confinement. 

The changes in the density profile as the kinetic energy 
becomes non-negligible can be discussed in terms of the 
dimensionless parameter 



ri = Na/\i_ 



(5) 



which is approximately the ratio between the TF-ground 
state energy per particle for 1=1 and the recoil energy 
Er = h^q^/2m (see Fig. 3a). When the product 77/ de- 
creases and becomes comparable to one, the radii of the 
condensate both increase with respect to those obtained 
in the TF-limit. This can be understood comparing the 
TF-energy, Etot ^ Nrj Er P , with the kinetic energy, 
Ekin ~ N En I, estimated using the TF-Gaussian solu- 
tions. 

We can compensate the kinetic energy pressure in the 
radial direction, when necessary, by adding external ra- 
dial confinement, e.g. by an appropriate choice of the 
focus of the laser beam. Fig. 3b shows the results of fix- 
ing the radial width Wr and decreasing the parameter 77. 
As soon as the kinetic energy associated with the density 
modulation, that is larger as N Er, exceeds the energy 



reduction caused by density modulation, e. g. 77/^ ^ 1, 
we find the oscillatory pattern to be washed out. For 77 
smaller than 2.8 at / = 1, the condensate is still self- 
confined along the longitudinal direction, with a similar 
value of the longitudinal size Az, but the density oscil- 
lations are almost washed out. Only for much smaller 
values of 7/, for instance ^ 0.02 at / = 1, does the kinetic 
energy also affect Az (long dashed line in Fig. 3b). 

Remarkably, when 7//^ ^ 1, the single condensate 
droplets overlap, self-organizing into a spatially- coherent 
long-range density modulation — a Bose "supersolid" |^ 
is formed. In this novel regime the atoms are coherently 
distributed among the wells allowing the establishment 
of a phase between the overlapping condensates. An ex- 
ample is given by the dot-dashed curve in Fig. 3b. 

Using the Gaussian ansatz in its time-dependent form 
we then obtain analytically in the small intensity 
limit the values for the radial and longitudinal oscilla- 
tion frequencies for a self-bound condensate in absence 
of the density modulation (/ <C rjP ^1): fi^ = 
2.88 {ER/hf {Pt]) and = 0.0148 {Er/H)^ {I*r]). It is 
noteworthy that, although the condensate radii are much 
larger than the laser wavelength (far-zone), the frequency 
fl^ of the compressional mode along the radial direction 
is close to the "plasma" frequency uj^ = ATrun{0)/m = 
9.86 [Er/Ti]'^ (I^v) introduced for the condensate in the 
near-zone "self-gravity" regimes [Q, n{0) being the peak 
density and u= (11/10) TiTray/q the "gravitational" cou- 
pling of the near zone laser-induced interatomic potential 
—u/r The longitudinal frequency Q,z is expected to 
be substantially modified in the presence of the oscilla- 
tory pattern. 

The spontaneous Rayleigh scattering at a rate Fray 
from the interaction-inducing laser beam leads to heat- 
ing and depletes the condensate. Treating the scattering 
of photons as an incoherent process, we expect the total 
energy of the condensate to increase as |l^ 
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The characteristic time in which the condensate may 
heat up is Thcat = {{d Etot / dt) / \Etot\)~^ ■ We estimate 
Thcat — -021^5 N {H/Er) I, using the relations for the 
small intensity limit obtained using the Gaussian ansatz 
{Etot/N = -0.71&9 tiErP), and compare it with the 
above expression for uip. The product TheatWp, that can 
be expressed in terms of the number of atoms per cubic 
wavelength as Thcatt^p = 0.11 (a/AL)i/2(Ain(0))3/2/-i/2^ 
proves that heating due to spontaneous Rayleigh scat- 
tering events does not provide a fundamental limitation 
on the observability of the discussed effects for sufficient 
density. 

For the density modulations discussed above, most of 
the photons are nearly elastically and coherently scat- 
tered. Their total cross section is given by N^fcm times 
the single-atom cross section, where the fraction fom is 
defined as an appropriate average over all the possible 
directions of scattering (in the Born approximation) 
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Here n(k) is the Fourier component of atomic density 
corresponding to k = 2q{sm{9)x , 0, cos{9)z), the momen- 
tum transfered by a single photon in the x-z plane and 9 
is the angle between the incident beam and the scattered 
direction. Only the (complementary) fraction (1 — fcm) 
of the rate of energy change (||), that is related to the in- 
coherent part of the scattering cross-section, is absorbed 
by the interatomic degree of freedom in the center-of- 
mass frame, thereby providing a partial suppression of 
heating. 

The center of mass is therefore subject to a con- 
stant radiation force that is enhanced by a factor N fcm 
by the density modulation. Correspondingly the elec- 
tromagnetic field is strongly coherently back-scattered 
(diffracted) . The radiation force can either shift the equi- 
librium position of the condensate, if it is located in a lon- 
gitudinal trap, or else accelerate it uniformly. In the lat- 
ter case the scattered light will be Doppler shifted. The 
effects described here are due to the same matter field 
interactions as those responsible for stimulated Rayleigh 
scattering [|l3| and collective atomic recoil (CARL) [Q. 
In fact, fcm will lead to similar effects. However, the 
change of atomic energy, unaccounted for by CARL, is 
the essence of our effects. 

A central prediction of this paper is that induced 
dipole-dipole forces result in a stationary density mod- 
ulation for the condensate ground-state. Density mod- 
ulations in a condensate can also occur in the presence 
of phonon excitations, as observed in the experiment by 
Inouye et al |^ that has demonstrated the superradiant 
Rayleigh scattering of a pulse of light by an atomic 
BEC. However, traveling phonons represent excitations 
of the system and so are quite distinct from the ground- 
state density modulation proposed here. The two sit- 
uations can, for instance, be distinguished by diffract- 
ing a nearly perpendicular probe laser with z-component 
of the wavevector, ~ 2q. Sound waves at frequency 
il — Vs2q, where Vs is the sound velocity, lead to a den- 
sity modulation of the form cos(2gz — fit) so that the 
first diffraction orders (i.e. Brillouin peaks) of the probe 
beam (wavenumber k), at angles 9 w ±2q/k, will be fre- 
quency shifted to ck ± fl By contrast, diffraction 
from the stationary density modulation discussed above 
will be elastic, with no frequency shift. 

An example of the experimental conditions required 
for the predicted effects involves the following parame- 
ters for a cloud of A'' ~ 10'^ sodium atoms: a circularly 
polarized laser beam, red-detuned by 1.7 GHz from the 
3Si/2 (F=l) -> 3P3/2 (F = 0,1,2) transition, gives the 
threshold for 1/r^ instability at « 525 mW/cm^ For 
such intensity, the acceleration of the center of mass is 
approximately « bOONfcm [m/s^], where fcm is a non- 
negligible fraction (for instance fcm ~ 0.1 for / 1 and 
fcm ~ 1 for / <C 1). It is possible to balance this force 
by combining longitudinal harmonic confinement with a 



magnetic field gradient, so that the condensate will not 
accelerate. 

To conclude, we have demonstrated a new regime 
of self-confined and self-organized ground-state density 
modulations in a BEC illuminated by a single circular 
polarized laser beam. We have shown that non-linear 
scattering of light may arise, even in the small-saturation 
limit, and have estimated the enhancement of the total 
"elastic" cross section. This regime is inherently possible 
even for a plane-wave laser, although it is facilitated by 
the radial focusing of the beam. 

This work was supported by the German-Israeli Foun- 
dation (GIF). 
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